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UNIVERSITY OF ILLINOIS 
DIGITAL CCMPOTER 

LIBRARY ROUTIN E LJ - 250 • 
By Dave Muller 

Automatic Linear Equation Solver With Programmed Checks 

and Calculation of Residues (SA90Z Only) 

Entire Jrogram 

Depends upon the condition of the equations to he solved* 

(a) 15 seconds to input program 

(b) «004k seconds to input coefficients, where k is the 
number of characters on the data tape 

(c) ,0017 n? seconds to solve equations. 

(d) »017(p+3)n seconds to punch results, where p is the 
number of decimal digits punched for each result. 

1^3 equations 

$iis program is an extension of the older routine 12 and is 

used in the same way as L2. It uses the drum whenever the 

number of ^equations exceeds 32, and is capable of solving 

as many as 1^3 equations. 

Read the program into the memory in the usual way. When the 

program tape has been read, place the data tape in the reader 

and restart; with the black switch* The results will be 

punched in a single column with a decimal point properly 

placed in each result. 

To solve the set of equations 

n-1 

j£> V X J +a in=° 

we proceed as follows t 

(a) Scale the coefficients by rows and also possibly by 
columns (usually simply by moving the decimal point) so 
that each coefficient is less than l/2. If the 
coefficients have been scaled by columns one must take 
account of this in interpreting the results* 

(b) Punch each scaled coefficient as a sign followed by 12 
or fewer decimal digits- The decimal point ia' assumed 
to follow the sigh so that +.016 and -*204 should be 
punched as +0l6 and -20^. 



•fi- 
fe) Terminate each rov a lQ , ..., a^ of coefficients 

by punching the character I. 
(d) Follow the last K by a sexadecimal character p which 
determines the number of decimal digits to be printed 
in the results* The character p can assume the values 
2, 3, . .., 9, K, S, ff. Spaces (5 holes) may be punched 
at will in the data tape. A leader of non*-sexadeclmal 
(5th hole punched) characters should precede the first 
coefficient. 
PROCLAIMED CHECKS If the program stops on the instruction FF 026 the drum has 

failed during triangularization. 

If the program stops on the instruction IT 027 the drum has 
failed during back substitution. Another attempt at back 
substitution may be made if this instruction is by-passed 
with the white switch. 

If the program stops on the instruction FF 028 the wrong 
number of coefficients has been read from the data tape • 
during the reading of the last equation. 
If the program stops on the instruction FF 029 the solution 
of more than 1^3 equations is being attempted. 
If the program stops on the instruction FF 02K before reading 
the coefficients the sum check on the master tape has failed, 
indicating that it was read incorrectly. 
CAIXTOATIOJf OF RESIDUES After completing the solution of a set of equations the 
^ computer will stop on the instruction 21*009. One may then 
read and solve another set of equations by throwing the 
black switch. Alternatively, after completing the solution 
of a set of equations one may calculate the residues by means 
" of a program for this purpose which is attached to the end 
of the L-7 tape. This program may be read by throwing the 
white switch rather than the black switch. The tape of 
coefficients should be reread using the black switch and 

the res idues 
n-1 



b ± « E 

will be computed and printed. 



i " J, "y x i + "in * * " ^ •••» ** 



-3- 



If the residues are not small compared to the roots one may 
conclude that the calculation of the roots has been carried 
out Inaccurately. On the other hand, the fact that the 
residues are small does not necessarily Indicate that the 
roots are accurate since the equations may be poorly conditioned 
or they may have been Ineptly scaled. 
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K 








00 3K 






3 


00 F 








00 3*2F 




r o 


k 


00 F 








00 9F 




Location of L 6 


5 


00 F 








00 22 TF 




Location of master routine 


6 


00 F 








00 171F 




Location of P 16 


7 


00 F 








00 270F 




Location of N12 


8 


00 F 








00 30S5 


• 






00 9K 






Modified Version of L 6 




1 1 

Word in Present Routine 


Word in L 6 


3 


26 S7 




26 S7 




26 S5 




L5 150L 


66 


32 63L 




32 63L 




kl F 




*U 1S9 


TO 


IA F 




IA 1S9 




kO F 




40 1S9 


115 


50 F 




50 S9 




26 87L 




26 87L 


120 


26 28S5 




50 lVfL 




7J Sk 




7J Sk 


129 


LO 1*9L 




LO 1*9L 




32 30S5 




36 SK 


138 


85 11F 




85 11F 




00 2560F 




00 SS 


139 


NO F 




NO 1S9 




00 F 




00 F 


1*3 


80 IF 




80 IF 




L5 F 




L5 1S9 

1 



L 7 



LOCATION 



1^5 

151' 

152 
161 







ORDER 



(P16) 
(N12) 



86 11F' 
00 2560F 

S2 67L 
L5 F 
NO F 
kO F 
75 F 
LO F 

00 171K 

00 270K 

00 227K 
L5 21S7 
k6 klL * 
10 20F 

k2 Jul 

Ll 3F 
56 7L 
L5 *UL 
LO 1^2S^ 
56 5L 
FF ^K)F 
LO 13^S4 
32 kL 
L5 150S4 
26 i+Sif 



NOTES 



PAGE 2 



86 11F 
00 SS 
S2 67L 

L5 2S9 
NO F 
to 2S9 
75 F 
LO 1S9 

Infraprint 
Infraput 



Plant r Q + n + 1 



Test for first row 



Test number of coefficients 



L 7 



LOCATION 


ORDER 


*OPES HUB 5 


-• 


7 


L5 4lL 
LO 136S4 








8 


LO 134S4 
40 146S4 








9 


LO 42L 
32 10L 








10 


FF 4lp 
L5 l46s4 








11 


L4 154S4 
L4 134s4 




•» 




12 


46 159S4 
42 159S4 








15 


L5 4lL 
40 142S4 








14 


46 135S4 
42 135S4 


• 






15 


LO 134S4 
42 39L 








16 


42 8?S4 
42 120S4 




Test for more than 144 coefficients 
and plant addresses 




17 


LO 134S4 
46 20S4 




• 




18 


L5 4lL 
L4 146S4 








19 


46 115S4 
L4 134s4 








20 


42 66S4 
46 70S4 • 








21 l 


42 70S4 • 
46 139S4 








22 


42 143S4 
42 161S4 








23 


L4 134s4 
42 133S4 








24 


42 151S4 










42 152S4 




_ . 





L7 



LOCATION 
25 

26 

27 
28 

29 
30 

31 

"32 

33 

34 

35 • 

36 

37 

38 

39 

40 

in 

42 



ORDER 

42 150S4 
L5.136S4 
*2 33L 
L5 40L 
46 34L 
26 6l 
L5 34L 
L4 40L 
46 34L 
50 147S4 
22 120S4 
81 7F 

00 17F 

L4 34l ' 

46 54L 

92 131F 

92 513F 

L5 F 

J4 F 

50 34L 

26 S6 
F5 33L 
42 33L 
LO 39L 

32 32L 
92 131F 
24 S4 
26 327F 
12 513F 
L5 F 
00 100F 
00 F 
LJ F 
40 F 
00 144F 
00 F 
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L 7 



Change scale factor 



Print results 



M 



Constants 



LOCATION 


ORDER 


NOTES mbr s 






00 309K 


1 


• 







50 500F 
50 L 




Read a row of coefficients 




1 


26 270F 
L5 17L 




-r 




2 


ho kL 
hi F 








3 


2L i+L 
S5 F 




• 




h 


50 F 
7^ 500F 




Compute a residue 




5 


Lh F 
kO F 








6 


L5 4l 
IA 143F 








7 


40 4L 
LO 14L 








8 


32 3L 


1 






92 131F 








9 


92 513F 
L5 F 








10 


5^ F 
50 10L 




f 




11 


26 171F 
F5 15L 




Print a residue 




12 


hO 15L 
LO 16L 








13 


36 L 










24 28L 









L7 
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Prepare parameters and print 



RESIDUES 



Print ROOTS and return to L-7 



LOCATION 
31 

52 

33 

35 



ORDER 

92 707F 
92 131F 
26 9F 
2^ 1?L 
L3 F 
3^ 28L 
FF ^2F 
26 28L 
SO 1408F 
06 2319F 
26 33L 
26 IN 
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L 7 



i 



